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ABSTRACT 

We investigate the growth over time of 20 massive (> 3 keV) clusters in a hydrody- 
namical simulation of the ACDM cosmology with radiative cooling. The clusters show 
a variety of formation histories: some accrete most of their mass in major mergers; 
others more gradually. During major mergers the long-term (temporally-smoothed) 
luminosity increases such that the cluster moves approximately along the Lx-Tx re- 
lation; between times it slowly decreases, tracking the drift of the Lx-Tx relation. We 
identify several different kinds of short-term luminosity and temperature fluctuations 
associated with major mergers including double-peaked mergers in which the global 
intracluster medium merges first (Lx and Tx increase together) and then the cluster 
cores merge (Lx increases and Tx decreases). At both luminosity peaks, clusters tend 
to appear spherical and relaxed, which may lead to biases in high-redshift, flux-limited 
samples. There is no simple relationship between scatter in the Lx-Tx relation and 
either recent or overall merger activity or cluster formation time. The scatter in the 
Lx-M and Tx-M relations is reduced if properties are measured within R 500 rather 
than i? v ir- 
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1 INTRODUCTION 

Clusters of galaxies are the largest virialized structures in 
the Universe and mergers between them are among the most 
energetic events. This paper presents results from resimula- 
tions of 20 massive clusters to investigate the effects of the 
merger history on the X-ray properties of the intracluster 
medium (ICM). 

In recent years high-quality observations of the ICM 
of distant clusters of galaxies have become available. This 
has enabled cosmological evolution studies to be carried out 
on clusters by lo oking at, for example, the high-redshift 
Lx-Tx relati o n (iFairlev et alJ l2000| [ Borgani et all 1200 it 
iHolden et all 120021 : IVikhlinin et alJ booa iNovicki et alJ 
2002). This paper will look in detail at the merger history of 
clusters of galaxies and how this affects their X-ray prop- 
erties. Previous o bservational studies of merging clusters 
of gala xies i nclude | M arkevitch et all i2002T) . iMauehan et al. 
<2003ft and iReiprich et alJ i2003h . and I Allen et all (l20o37 
found evidence that the hottest cluster known (RX J1347.5- 
1145) is undergoing a merger. We find examples in our sim- 
ulated clusters that mimic these observed clusters. 

The merger rates of simulated Cold Dark Matter 
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(CDM ) h aloes were first inve stigated by lLacev fc Cold 
l|l994ft and lNavarro et all dl995Ti . The results were found to 
be in a greement with the ana lytic al model oflL acey fc Cold 
ll 1993ft . I Wechsler et alJ J2002h and lZhao et alJ (12003ft looked 
at the effects of mergers and accretion on the structure of 
individual CDM haloes and found that halo concentration 
tends to increase with increasing formation redshift. 

To investigate the effect s of m ergers on the ICM, 
iPearce. Thomas fc Couchmanl jl994ft simulated the col- 
lision of pairs of relaxed cluster haloes and com- 
pared these with earlier dark-ma tter only simulations 
iPearce. Thomas fc Couchmanlll993ft . They found that the 
entropy structure of both the dark matter and the gas is 
relatively unchanged during the merger, although there is 
a small tendency to transfer energy from the former to the 
latter in the core of the system. The time evolution of the 
X-ray properties o f merg ing clusters were investigated by 
iRitchie fc Thomas! i2002ft in higher resolution simulations 
that also included radiative cooling. They showed that the 
X-ray temperature and luminosity can temporarily increase 
by a large factor during a merger; the central entropy in- 
crease after the merger was found to be a strong function of 
subcluster mass and impact parameter. 

The simulations described above started with isolated, 
relaxed clusters and investigated different impact parame- 
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Table 1. Simulation parameters for the three runs: name; ef- 
fective number of particles; particle mass in units of h~ 1 MQ\ 
softening in units of h~ 1 kpc and particle species simulated. 

Name Particles Mass Soft. Particle types 

LDM 160 3 1.9 X 10 11 50 Dark matter 

HDM 320 3 2.4 X 10 10 25 Dark matter 

HGAS 2 X 320 3 2.4 X 10 10 25 Dark matter + gas 

ters and mass ratios. However, real clusters in a cosmologi- 
cal environment are more complicated. They will have some 
amount of substructure, may be rotating and could collide 
with more than one subclump at a time. Therefore to find 
out what kind of mergers clusters tend to undergo and the 
effect this has on the complex ICM, full cosmological simu- 
lations need to be carried out. A first step in th is direction 
was undertaken bv lEke. Navarro fe Frenkl jl99ct) . They un- 
dertook resimulations of individual clusters extracted from 
a dark-matter simulation and showed that the evolution of 
bulk properties of clusters varied from cluster to cluster. 
The simulations that we describe in this paper are similar 
in spirit but have more particles, include radiative cooling 
of the gas component, and have a much higher time reso- 
lution with which to follow the development of the X-ray 
properties. 

We should also men tion a couple of other recent stud- 
ies. iRandall et alJ J2002) looked at the observational bias 
that can be induced by the temporary enhancements in lu- 
minosity and temperature in merging clusters on determina- 
tions of as and th e number density of high-redshift clusters. 
iMotl et alJ i2003l) undertook radiative simulations of a sam- 
ple of two clusters and showed that the traditional model of 
smooth accretion onto clusters is inaccurate in that clusters 
accrete gas in subclumps which bring precooled gas direct 
to the core of a cluster. Our work supports both of these 
ideas. 

The rest of this paper is organized as follows. In Sec- 
tion 2 the properties of the simulations and the resimulation 
technique will be described. Section 3 will investigate the 
mass, luminosity and temperature evolution of the resimu- 
lated clusters, with particular types of event being identified 
and examples of these looked at in detail. The link between 
scatter in the scaling relations and merger history or sub- 
structure will be investigated in Section 4. Finally, in Sec- 
tion 5, we summarise our results and compare the features 
seen in the simulations to observations. 



2 THE SIMULATIONS 

In a previous paper, iMuanwong et all l)2002l hereafter 
MTKP02), X-ray scaling relations were presented for sim- 
ulated catalogues containing over 400 groups and clusters. 
However, because of the volume simulated, only 7 of these 
had virial temperatures in excess of 3 keV. This study sim- 
ulates a further 20 such clusters by resimulating regions 
within a larger box, with the high-resolution regions having 
simulation parameters identical to that of the earlier run. 

Three simulations were undertaken with parameters as 
shown in Table The effective number of particles is the 



number of particles needed to fill the entire simulation vol- 
ume at the highest resolution in the box. In the simulation 
HGAS there were equal numbers of gas and dark matter 
particles in the high-resolution regions. The mass resolution 
for HGAS refers to the sum of the two species. 



2.1 The resimulation technique 

We first generated a low-resolution run (LDM) with 160 3 
dark-matter particles. This was run to completion (from 
z — 49 to z = 0) and the 20 largest clusters identified using 
the technique described in MTKP02. All the particles within 
2 virial radii of the cluster centres were identified. The lo- 
cations of all these particles on the initial comoving grid 
(i.e. before the Zel'dovich approximation is applied) were 
noted and, along with the neighbouring grid points, were 
flagged as requiring high resolution. Just under 4 per cent 
of the box was included in the high-resolution mask. 

A second set of initial conditions was then generated on 
a 320 3 grid. The original waves used in the LDM simulation 
were used again with additional new waves generated for 
frequencies between the LDM & HDM Nyquist frequencies. 
The high- and low-resolution initial conditions were then 
combined with the former being used within the previously- 
identified mask and the latter elsewhere. This resulted in 
5 191 535 dark-matter particles (as compared to 4 096 000 in 
the original low- resolution simulation). 

Finally a third set of initial conditions was gener- 
ated that included both gas and dark matter in the high- 
resolution regions. To do this, a gas particle was added at the 
location of each high-resolution dark-matter particle. The 
masses of the two were adjusted to give the same total mass 
as before and to give the universal dark matter to baryonic 
matter mass ratio. The dark-matter and gas particles were 
given identical displacements and velocities. The total num- 
ber of particles for this run was 6 443 575. 

The hydrodynamics simulation used parameters identi- 
cal to that for the Radiative model of MTKP02. In partic- 
ular, it included radiative cooling with a time-varying met- 
alicity Z(t) = 0.3 (t/to) Zq, where to is the current age of 
the universe (approxima tely 12.8 Gyr for this cosmo logy), 
using cooling tables from ISutherland fc Dopital lll993T) . Nei- 
ther preheating nor feedback were included in this run. 

The simulations were run on 64 processors on the 
Cray T3E at the EPCC, using a parallel form of the 
AP 3 M SPH code HYDRA dCouchman et al.l 11995b 1 . The 
box was 200 ft -1 Mpc across and used cosmological param- 
eters: h = 0.71; a 8 = 0.9; (l m = 0.35; D. v = 0.65 and 
Ob = 0.038. The shape parameter for the fluctuation spec- 
trum was fixed at r = .21 using the fitting function of 
iBond fe Efstathioul <ll984Tl . The gravitational softening was 
set to 100/(l+z) /i _1 kpc until z = 3, after which it was held 
fixed at 25/i~ 1 kpc. The initial redshift was z = 49 and the 
highest-resolution run required 3706 timesteps to evolve to 
z = 0. 



1 HYDRA is available from http://hydra.susx.ac.uk/ 
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Table 2. Cluster properties: identifying number; virial radius 
in units of /i — 1 Mpc; virial mass in units of 10 14 ft, _1 MQ; virial 
temperature in units of keV/fe; sub-structure statistic; soft- 
band, emission-weighted temperature excluding emission within 
50/i kpc of the cluster centre in units of keV/fc; bolomctric 
luminosity estimated from soft-band emission, excluding emis- 
sion from within 50/i -1 kpc of the cluster centre, in units 
10 44 /i — 2 ergs — 1 \ formation expansion factor as fit to the M v - lT 
evolution-see Section 3.1. The clusters marked in bold are ones 
which will be examined in detail below. 



id 


R ■ 


M v i r 


j" . 


T x 


L x 


Sub. 




1 


2.41 


17.85 


9.62 


7.15 


9.22 


0.11 


0.671 


2 


1.80 


7.49 


4.99 


4.68 


1.37 


0.05 


0.529 


3 


1.79 


7.33 


4.63 


3.41 


1.13 


0.14 


0.718 


1 


1.78 


7.26 


5.09 


5.09 


1.73 


0.02 


0.688 


5 


1.67 


6.02 


3.75 


3.07 


0.43 


0.12 


0.890 


6 


1.62 


5.46 


4.04 


4.09 


1.04 


0.07 


0.557 


7 


1.62 


5.44 


3.61 


3.10 


0.38 


0.19 


0.756 


8 


1.61 


5.39 


4.27 


5.24 


2.63 


0.03 


0.443 


9 


1.61 


5.34 


4.70 


3.94 


1.93 


0.15 


0.845 


10 


1.60 


5.25 


4.09 


4.33 


1.36 


0.02 


0.459 


11 


1.60 


5.22 


4.15 


4.19 


1.74 


0.08 


0.731 


12 


1.59 


5.13 


4.37 


4.58 


2.19 


0.03 


0.565 


13 


1.58 


5.06 


5.38 


6.18 


3.12 


0.09 


0.487 


14 


1.56 


4.83 


3.83 


3.77 


0.52 


0.25 


0.774 


15 


1.54 


4.71 


3.41 


4.11 


1.27 


0.06 


0.535 


16 


1.54 


4.71 


3.15 


3.00 


1.14 


0.11 


0.777 


17 


1.53 


4.61 


4.59 


4.19 


1.92 


0.06 


0.501 


18 


1.53 


4.58 


3.62 


4.27 


0.84 


0.03 


0.699 


19 


1.51 


4.43 


3.48 


4.38 


0.99 


0.05 


0.480 


20 


1.35 


3.17 


3.13 


3.60 


1.02 


0.06 


0.401 



2.2 Cluster properties 

The properties of the resimulated clusters at z = are 
listed in Table [5] The virial radius is defined as the radius 
of a sphere, centred on the densest dark-matter particle, 
that encloses a mean density of 317 (specifically 178 r2m°' 55 > 
lEke et all 1998) times the background density (111 times 
the critical density). The virial mass, M v ir, is the mass en- 
closed by this sphere and the virial temperature is the mean 
specific energy (kinetic plus thermal) of the dark matter 
and gas, multiplied by fimn/k, where fimn = 10 _24 g is 
the mean molecular mass and k is the Boltzmann constant. 
The X-ray temperature is weighted by emission in the soft 
X-ray band (0.3-1.5 keV), excluding emission from within 
50/i _1 kpc (i.e. physical, not comoving) of the cluster cen- 
tre (hereafter referred to as cooling- flow corrected emission). 
The soft-band X-ray luminosity is converted into an esti- 
mated bolometric luminosity using the procedure described 
in MTKP02. The substructure statistic is denned as the sep- 
aration between the centroid of the mass and the location of 
the densest dark-mat ter particle, in units of the virial radius 
llThomas et al Jl998h . The formation expansion factor, <jf, is 
a characteristic formation epoch for the cluster as denned in 
Section 3.1. 

Note that, after resimulation, the virial mass and tem- 
perature of the final cluster in the list dropped significantly. 
This is because the cluster is in the process of merging with 
a subcluster that lay just inside the virial radius in the orig- 
inal low-resolution run but has moved just outside it in the 



high-resolution run. Even so, this cluster still has over 13 000 
particles each of gas and dark matter within the virial ra- 
dius at the final time. The largest cluster has approximately 
73 000 particles of each type at the end. 

2.3 Testing 

The mass of dark-matter particles in the low-resolution re- 
gions exceeds that of those in the high-resolution regions by 
a factor of ~9 (and of the gas particles by a factor of ~74). 
It is therefore necessary to ensure that the clusters are not 
affected by these high mass particles which could cause un- 
physical two-body relaxation (since the softening is the same 
for all the particles). The clusters were examined to discern 
if there were any low-resolution particles within their virial 
radii at any output time. It was found that only one of the 
clusters presented was at all affected. This was Cluster 13 
which temporarily has up to 2 low-resolution dark-matter 
particles within the virial radius during a late-time merger. 
We do not expect this to affect the results and we choose to 
ignore it. 

As a test of our method the clusters were compared 
to those from the Radiative simulation of MTKP02. A 
comparison between the temperature-mass relation and the 
luminosity-temperature relation for the two is shown in 
Fig. □ To be consistent with MTKP02, the extent of the 
clusters is defined by a sphere that encloses a mean den- 
sity of 200 times the critical density. It is clear that the new 
clusters are consistent with the previous relations but extend 
them to higher mass and temperature. 

One surprising result, however, is that the number of 
clusters above a virial temperature of 3 keV in the new runs 
is far fewer than 8 times the number in the original MTKP02 
simulation (20 as compared to 7) which was one-eighth of 
the volume. To test our normalization, we compare in Fig. 
the mass -function of the simu lated clusters with that pre- 
dicted bv I Jenkins et alJ i200ll) from a compilation of a large 
number of A r -body simulations. The 20 high-resolution clus- 
ters are contained in the final three bins of the plot (except 
for the most massive cluster that lies off the right-hand-side 
of the plot). In the same mass-range, 3 clusters are expected 
in the MTKP02 simulation but 6 are found. We put this 
down partly to chance and partly due to the fact that simu- 
lations do not sample waves correctly on scales comparable 
to the box-size. 



3 THE TIME EVOLUTION OF MASS, X-RAY 
LUMINOSITY AND TEMPERATURE 

In this section we look at the change in the mass, X-ray 
luminosity and X-ray temperature of the clusters as they 
evolve. For simplicity and for ease of physical interpretation, 
we use the bolometric (i.e. not the soft-band used above) 
emission from within the virial radius in the rest-frame of 
the cluster. 

The time resolution is defined by the light crossing time 
of half the box, that is 0.5 x 200/i~ 1 Mpc x a/c « 4.6 x 
10 8 x ayr, where c is the speed of light and a = 1/(1 + z) 
is the expansion factor. These times were chosen to match 
previous simulations which were used to obtain light-cones. 
This gives sufficient resolution to crudely resolve the mergers 
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Figure 1. Cooling-flow corrected Tx — M200 (upper panel) and 
Lx — Tx (lower panel) relation in the soft X-ray band (luminosity 
converted into an estimation of the bolometric luminosity) for 
the clusters from MTKP02 (squares) and those from simulation 
HGAS (stars). The lines have a slope of 0.59 (upper panel) and 
2.80 (lower panel). 



for the clusters as a whole, although finer time-resolution 
would have allowed us to examine accretion into the cores 
of the clusters in more detail. 



3.1 The growth of mass 

Fig. shows the time-development of the mass enclosed 
within -Rvir- The most-massive cluster is located in the top- 
left panel of this figure (and also in Figs 2] 03 & IHJ , with 
successively smaller clusters reading from left-to-right and 
then top-to-bottom as numbered. The plots go back to a 
time (a = 0.35) when the smallest cluster has a virial mass 
of 1.5 x 10 13 /i _1 Mq, corresponding to 600 particles of each 
species. Also shown is M500, the mass enclosed by a sphere 
within which the mean density is 500 times p c o(l + z) 3 where 
p C Q is the critical density at z = 0. 

It is apparent from this figure that there are times when 
the mass of a cluster increases smoothly and other times 
when it undergoes sudden jumps; we designate these as mi- 
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Figure 2. The mass function of clusters in the LDM simulation 
(dashed line) and in the simul ations of M TKP02 (dotted line) 
compared to the prediction from Jenkins et al. 12001) (solid line). 
The most-massive cluster lies off the right-hand-side of the plot 
and has not been included in the figure as there are numerous 
empty bins between it and the second most-massive cluster. 



nor and major mergers, respectively. Although the division 
between the two is somewhat arbitrary, it is a useful one, 
as we will show below that X-ray cluster luminosity evolves 
quite differently in the two regimes. For the purposes of this 
paper we define the onset of a major merger as occurring 
when the jump in log 10 M v i r between one output time and 
the next is 0.08 or more (i.e. a factor of 1.2). Visual inspec- 
tion of the luminosity and temperature evolution suggests 
that the effects of a major merger last for three successive 
output times, approximately 1.4 x 10 9 x ayr. This roughly 
corresponds to the dynamical time of the clusters, although 
it does not scale in quite the same manner (oc a instead of 
oc a 3//2 ). The jump criterion of 0.08 is chosen such that, aver- 
aged over all the clusters, about half of the mass is accreted 
during major and half during minor mergers. On average, 
about 25-30 per cent of the cluster outputs correspond to 
periods of major merger activity. 

From Fig. it can be seen that clusters exhibit a wide 
variety of formation histories. Some (e.g. clusters 9 or 16) 
aquire most of their mass in major mergers, while others 
(clusters 2 and 20) undergo no major mergers at all. The 
percentage increase of logarithmic mass during major merg- 
ers is show in the figure and summarized as a histogram in 
Table 01 

It is common in the literature jLacev fc Cold Il994l : 
iNavarro et al.l ll997') to tacitly assume that cluster proper- 
ties are fixed at the time of formation and remain unchanged 
since then. However, only in a minority of cases is it possi- 
ble to assign a particular 'formation time', associated with 
a major merger, to the cluster. An altern ative approach has 
been suggested bv lWechsler et all (120021) who fit models of 
the form M(a) = Mo e 2af(1 ~ to the data. The point where 
the slope of the relation equals 2 defines a characteristic for- 
mation epoch for each cluster, as listed in Table |5] We shall 
investigate later, in Section 4, whether the X-ray properties 
of clusters are correlated with their formation time. 

The growth in M V [ T is relatively monotonic although 
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Figure 3. M v ; r (thick line) and M500 (thin line) versus expansion factor for the clusters. The dashed curves are a one-parameter fit 
to the functional form M(a) = Mo e 2at a' introduced bv lWechsler et all i2002t) . where af is listed in Table 2. The asterisks indicate 
events which will be investigated in detail below. The lines along the bottom of each box indicate the times during which the cluster is 
deemed to be undergoing a major merger (as defined in the text). The numbers in parentheses indicate the fractional logarithmic mass 
accreted during major mergers. 
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Table 3. Number of clusters accreting various fractions of their 
material through major mergers. Also show is the mean angle 
(with respect to positive Tx evolution and no Lx evolution) in 
log space which the clusters move accross the Lx — 1% plane (see 
Fig.®. 

Fraction accreted >0% >25% >50% >75% 

during mergers <25% <50% <75% <100% 

Number of clusters 4 4 8 4 



mean Lx-Tx relation as gradually drifting to lower normal- 
ization as time increases. Major mergers do not affect this 
drift but shift individual clusters along the relation to higher 
temperature and luminosity. 

Simple scaling arguments suggest that L oc MpT 1 ' 2 
where M is the cluster mass, p the density and T the tem- 
perature. (Note that the dominant contribution to the flux 
comes from the region where p oc r _3//2 so this defines the 
characteristic size to use in the scaling relation.) For the 
population as a whole one would expect the luminosity to 
decrease with time, at fixed mass, as the characteristic den- 
sity decreases roughtly in proportion to the mean density of 
the Universe. However, for individual clusters, the increase 
in mass more than makes up for this. However, as we have 
seen above, the increase in luminosity is not smooth. Dur- 
ing major mergers the mass increases abruptly. There will 
be an influx of cool, low-entropy gas and much of this will 
reach the core, causing a long-term increase in luminosity. 
Between major mergers the mass continues to increase but 
the influx of low-entropy gas is insufficient to replenish gas 
lost by radiative cooling and the core density and luminosity 
decrease. This will be looked at in more detail in a future 
paper. 



Mean angle of Lx-Tx drift -47° 14° 21° 38° 
Mean dlog L x /dlog T x -1.07 0.25 0.38 0.78 

there is sometimes a slight decline after a merger. The rises 
in M500 occur later and are often followed by a decline as 
the central regions of a cluster settle down after the merger. 
For this reason, it may be thought that M v i r is a more useful 
measure of the mass of a cluster. However, it will be shown 
in Section 4 that the X-ray properties are more closely cor- 
related with M500, because the bulk of the X-ray emission 
originates in this central region. 



3.2 Long-term trends in X-ray properties 

Figs 2] & HI show the time-development of the bolometric 
luminosity and emission-weighted temperature of the clus- 
ters. The thick line shows the total emission and the thin 
line shows cooling-flow corrected emission. 

These figures show that each cluster's X-ray tempera- 
ture and luminosity are undergoing continual fluctuations 
associated with both major and minor mergers. These tend 
to obscure the long-term trends, and so we have defined 
smoothed versions of these plots in which the average pro- 
files are defined as the median from each set of 5 successive 
output times. 

Fig.|Slshows the trajectory of the smoothed X-ray prop- 
erties for each cluster in the luminosity-temperature plane. 
Also displayed are the fractional logarithmic mass accreted 
during major mergers, and vectors indicating the overall 
drift of the cluster across the Lx — Tx plane. It can be seen 
that the clusters that accrete only a small amount of mass 
during major mergers tend to move down and to the right 
(i.e. they get less luminous with time), whereas those that 
accrete most of their mass in this way move up and to the 
right (i.e. they get more luminous). This is confirmed by the 
statistics in Table |!|] which show the average direction of mo- 
tion across the Lx-Tx plane in four bins corresponding to 
different fractional logarithmic mass accreted during major 
mergers. 

For the majority of their lives, between major mergers, 
clusters grow dimmer and heat up (on average increasing 
their temperature by 0.20 in the log with a gradient in the 
log Lx-log Tx plane of -1.46) so that the normalisation of the 
Lx-Tx relation decreases with time. During major mergers, 
clusters get brighter and hotter (on average increasing their 
temperature by 0.26 in the log with a gradient of 1.54). If 
we allow for the time the mergers take then this motion 
is roughly parallel to the mean relation (L x oc T x 2 8 ). This 
would suggest the merger history of clusters does not induce 
significant scatter in the relation. Thus we can think of the 



3.3 Short-term fluctuations in the X-ray 
properties 

We next turn from the long-term behaviour of the smoothed 
temperature and luminosity profiles to look at the fluctua- 
tions during mergers. These exhibit a variety of behaviours 
that we shall illustrate with specific examples drawn from 
our simulations. Figs [7| M ED & HI show maps of several 
example clusters at different expansion factors as indicated 
in the ca ptions. The maps are generated in the same man- 
ner as in lOnuora. Kav fc Thomas! l)2003T) . The half- width of 
each map is equal to the virial radius and the emission is 
projected along the line-of-sight on each side of the clusters 
to a depth equal to twice the virial radius. The contours 
show X-ray surface brightness and are separated by 0.2 dex; 
the colours indicate emission-weighted temperature; and the 
arrows show the mass- weighted bulk flow, normalized to the 
highest velocity grid point. The corresponding evolution in 
the luminosity-temperature plane is shown in Figs ISl tTUl 1121 
& I14I in order to show the relative size of fluctuations, these 
figures all have the same dynamic range. 



3.3.1 Slow accretion 

Even when not undergoing major mergers, clusters still 
show fluctuations in their luminosity and temperature evo- 
lution. However, these are reduced if emission from within 
50 ft _1 kpc of the cluster core is excluded (see Table^J. What 
is happening is that the clusters are growing by accretion of 
small subclumps and the fluctuations are due to emission 
from low-entropy gas that makes its way to the cluster core 
and cools rapidly. This sometimes causes a temporary de- 
crease in the temperature of the cluster core, such as the two 
dips in the temperature of Cluster 10 at expansion factors 
of 0.60 and 0.67. 

Cluster 2, the second most-massive cluster in our sam- 
ple, has accreted all its mass during minor mergers. Fig. 
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Figure 4. Bolomctric luminosity versus expansion factor for the clusters. The thick line shows the total emission; the thin line shows 
cooling-flow corrected emission. The asterisks indicate events which will be investigated in detail below. The lines along the bottom of 
each box indicate the times during which the cluster is deemed to be undergoing a major merger. 
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Figure 5. Bolometric emission-weighted temperature versus expansion factor for the clusters. The thick line shows temperature calculated 
using the total emission; the thin line using cooling-flow corrected emission. The asterisks indicate events which will be investigated in 
detail below. The lines along the bottom of each box indicate the times during which the cluster is deemed to be undergoing a major 
merger. 



© 2003 RAS, MNRAS 000.IT1IT71 



The merger history of clusters 9 
r 




0.5 1 0.5 
lo 9io (kTx_ bol /keV) 

Figure 6. The time-development of the smoothed bolometric luminosity versus emission-weighted temperature for the clusters. The 
arrows show the overall drift of the cluster in the Lx — Tx plane. The bracketed numbers indicate the fractional logarithmic mass accreted 
during major mergers. 
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-0,4 -0.2 0,0 0,2 0.4 0,6 0.3 1,0 
log CkT/keV) 

Figure 7. Surface-brightness and emission-weighted temperature map for Cluster 2, as described in the text. The panels are at a = 0.544, 
0.560 & 0.577 so there is about 2.6x 10 8 yr between panels. The virial radius increases from 0.84 h~ 1 Mpc at the first panel to 0.88 /i -1 Mpc 
at the third panel. 



Table 4. Root-mean-square change in the log of the bolomctric 
luminosity and temperature between one output time and the 
next, excluding times during major mergers. 

Property Total Cooling-flow corrected 

Luminosity 0.120 0.090 
Temperature 0.069 0.059 




kT x /keV 

Figure 8. Evolution in the Lx-7x plane of Cluster 2 over the 
time shown in Fig. |7| The panels in Fig. |7| correspond to the 
numbers 1—3 here. The thick line plots properties for the whole 
cluster and the thin line plots cooling-flow corrected Lx an d Tx- 
The straight solid line shows a power-law relation of the form 
LxccTx 28 . 



shows snapshots at a fa 0.56 in which the infall of small sub- 
clumps is clearly visible. The substructure statistic is less 
than 0.1 throughout the evolution of the cluster, indicating 
that the subclumps are too small to significantly displace 
the centre of mass. Even so, the accretion gives rise to small- 
scale scatter as seen in Figs ^ [K] & |H] At any particular time, 



most clusters show this kind of lo w-level activity sim ilar to 
the minor mergers investigated bv lMotl et al] J2003). 

3.3.2 Single-peaked major mergers 

The most-common behaviour seen in major mergers is a 
temporary upward fluctuation in luminosity that is usually 
accompanied by an increase in temperature (the residuals of 
Lx and Tx have a Spearman rank coefficient of 0.48 which 
corresponds to a correlation with a significance of over 99.9 
per cent). Many examples of this can be detected in Figs|l|& 
|^| These fluctuations are not restricted to the core and cor- 
respond to more violent mergers that boost the luminosity 
of the whole cluster through compression of the intracluster 
medium. 

Figs. |5] & 1101 show an example from Cluster 8 at 
a ~ 0.63. This is caused by an approximately equal-mass 
merger (of clusters with masses 1.82 x 10 14 ft" 1 M Q and 
1.41 x 10 14 /i" 1 Mg) that produces a planar compression 
front perpendicular to the direction of the merger. The pla- 
nar nature of the front can be more or less obvious than 
that shown here, depending upon the impact parameter of 
the collision, the degree of substructure and the viewing an- 
gle. 

3.3.3 Double-peaked major mergers 

Sometimes a merger will be associated with a double peak in 
the luminosity, the first associated with an increase in tem- 
perature and the second a decrease. This happens when the 
cores of the main cluster and the in-falling subcluster do not 
merge directly but orbit about each other for a while before 
doing so. Unfortunately the time resolution of our outputs 
is insufficient to determine whether this is the generic be- 
haviour. 

Fig. II II shows a major merger of this type in Cluster 13 
at a ps 0.6. In the first two panels the subclump approaches 
and in the third panel the clusters collide and the temper- 
ature and luminosity peaks. At the time of maximum lu- 
minosity, however, the surface brightness contours are rel- 
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Figure 9. Surface-brightness and emission-weighted temperature map for Cluster 8, as described in the text. The panels are at a = 0.595, 
0.631 & 0.670 so there is about 5.8x 10 8 yr between panels. The virial radius increases from 0.68 h~ 1 Mpc at the first panel to 0.76 h~ 1 Mpc 
at the third panel. 
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Figure 11. Surface-brightness and emission-weighted temperature map for Cluster 13, as described in the text. The first panel is at 
a = 0.560 and the last is at a = 0.650 so there is about 2.8 X 10 s yr between each panel. The virial radius increases from 0.72 h' 1 Mpc 
at the first panel to 0.95 h _1 Mpc at the sixth panel. 



atively round and there is little evidence from the surface- 
brightness alone that a merger is taking place (the asym- 
metry that is visible in the figure is from in-falling matter 
that does not contribute significantly to the luminosity and 
temperature variation). The cores have not lost all relative 



momentum however and so they reseparate and do not fully 
merge until the final panel. This brings about a second peak 
in the luminosity which is only observed in the core (see 
Fig. ll2|l . Since this is relatively cool, the X-ray temperature 
drops at this point before stabilizing at a higher value. Note 
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Figure 10. Evolution in the Lx-Tx plane of Cluster 8 when 
it undergoes the merger shown in Fig. [5] The panels in Fig. [5] 
correspond to the numbers 1-3 here. The circles correspond to 
times when the cluster has a substructure statistic greater than 
0.1; other times are marked with triangles. The thick line plots 
properties for the whole cluster and the thin line plots cooling-flow 
corrected Lx and Tx. The straight solid line shows a power-law 
relation of the form Lx oc Tx 2 ' 8 . 




kT x /keV 

Figure 12. Evolution in the Lx-Tx plane of Cluster 13 when 
it undergoes the merger shown in Fig. 1111 The panels in Fig. 1131 
correspond to the numbers 1-6 here. The circles correspond to 
times when the cluster has a substructure statistic greater than 
0.1; other times are marked with triangles. The thick line plots 
properties for the whole cluster and the thin line plots cooling-flow 
corrected Lx and Tx. The straight solid line shows a power-law 
relation of the form Lx oc Tx 2 ' 8 . 

that, at the point the cores merge, the cluster again looks 
fairly spherical and relaxed. 

This is the behavio ur observed in the toy models of 
iRitchie fc Thomas! 1)20021) . The largest peak in the luminos- 
ity in these simulations was when the clusters first collided. 
Then the luminosity dropped as the core expanded again 
and then increased to a stable level as the core finally set- 
tled. This behaviour was more pronounced in off-axis colli- 
sions where the cores missed each other on the first pass. 




2 5 
kT x /keV 

Figure 14. Evolution in the Lx-Tx plane of Cluster 11 when it 
undergoes the merger shown in Fig. 1131 The panels in Fig. ll3l cor- 

respond to the numbers 1-3 here. The thick line plots properties 
for the whole cluster and the thin line plots cooling-flow corrected 
Lx and Tx- The straight solid line shows a power-law relation of 
the form L x oc T x 2 ' 8 . 

3.3.4 Major mergers leading to a permanent increase in 
luminosity 

Finally, there are major mergers that cause a permanent 
jump in luminosity and temperature after the initial fluc- 
tuation associated with the merger has died away. Half of 
the mergers involve a jump in the smoothed Lx (i.e. a per- 
manent jump in Lx) of over a factor of 2 and half of the 
mergers involve a jump in the smoothed Tx of more than a 
factor of 1.4. 

Fig. 1131 shows images before, during and after one such 
merger which Cluster 11 undergoes at a « 0.6. There is 
a small amount of hot, compressed gas as the subclump 
approaches but not so obviously planar as in some other 
mergers. At the point of merging the cluster looks fairly 
relaxed in this projection although the substructure statistic 
is greater than 0.1 throughout. 

Fig. 1141 shows that Lx and Tx increase together ap- 
proximately parallel to the Lx-Tx relation in an elongated, 
clockwise ellipse. This is because the luminosity increases 
when the clump crosses 7? v ir, before it has had a chance to 
interact strongly with the ICM and raise the temperature 
of the gas. A similar effect is seen in Fig. I1UI and the first 
(main) merger in Fig. 1121 



4 SCATTER IN THE X-RAY SCALING 
RELATIONS 

To a greater or lesser degree all of the X-ray scaling relations 
show scatter both in observations and simulations. Some of 
this will be due to observational errors and resolution effects. 
However, as clusters have differing merger histories and are 
not self-similar, it is apparent that some, perhaps most, of 
this scatter is physical. As Figs |3] & show, much of the 
scatter in luminosity and temperature comes from gas within 
50 h~ kpc of the cluster centre. For this reason, we omit the 
central gas when plotting the scaling relations. 
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Figure 13. Surface-brightness and emission-weighted temperature map for Cluster 11, as described in the text. The panels are at 
a = 0.577, 0.613 & 0.631 so there is about 5.5 X 10 8 yr between the first 2 panels and about 2.9 X 10 s yr between the second and third 
panels. The virial radius increases from 0.89 h^ 1 Mpc at the first panel to 1.04 h~ 1 Mpc at the third panel. 



Fig. 1151 shows the current-day X-ray temperature-mass 
relation within 7? v ir (top panel) and Rsoo (bottom panel). 
The straight line is the best-fit relation to a straight line of 
slope 0.60 in log space. This slope was chosen to agree with 
that from the clusters in MTKP02 which cover a much wider 
dynamic range in mass. We have investigated a wide variety 
of properties to try to find the main cause of the scatter. 
These include the following: 

• The most obvious source of scatter is the fluctuations 
in temperature associated with mergers, i.e. the difference 
between the actual temperature and the smoothed temper- 
ature. However, at the present day this is of a magnitude 
(approximately 0.03 dex) that is too small to have much of 
an effect. We have looked for a correlation in these fluctu- 
ations and deviations from the temperature-mass relation, 
and find none. 

• The fractional logarithmic mass accreted during major 
mer gers is also uncorrelat ed to scatter. 

• IWechsler et alJ (2002) found a correlation between clus- 
ter formation time and the concentration of the dark mat- 
ter. We used their method to assign characteristic formation 
epochs to each cluster (see Table |5J but again found no cor- 
relation with scatter. 

• The rate of increase of mass of both the cluster as a 
whole and of the core mass, averaged over different time 
periods before the present. This was also unsuccessful. 

In the end, only one of the statistics that we looked 
at gave a strong correlation with scatter from the X-ray 
temperature-mass relation and that was the degree of sub- 
structure. The substructure statistic, defined as the sepa- 
ration between the position of the centroid and the dark- 
matter density maximum in units of the cluster radius, is 
listed in Table [5] (Note that this statistic is a function of 
the cluster radius; the values in the table are for properties 
averaged within the virial radius.) The substructure statistic 
has several advantages over some of the other measures that 
we tried: it is relatively simple to calculate, it depends only 
upon the properties of the cluster at the present day, and it 
is closely related to some observable quantity (for example 



Table 5. Scatter about the fits to the plots in Figs 1151 to 1171 
This is defined as the root-mean-square deviation in the vertical 
direction in dex after allowing for the one free parameter in the 
fit. 

Relation Scatter in R v i T values Scatter in H500 values 

T x - M 0.088 0.057 

L x — M 0.25 0.19 

L x -T x 0.16 0.15 



ISchuecker et alJ l)200rJl observed substructure in 52 ± 7 per 
cent of clusters from the REFLEX+BCG sample) . 

In Fig. 1151 clusters with a substructure statistic greater 
than 0.1 are plotted as circles. It is immediately apparent 
from the upper panel that a large part of the scatter is 
related to substructure (formally, the distance below the 
line correlates to substructure with a Spearman rank co- 
efficient of 0.485, a significance of more than 95 per cent). 
It is interesting that most clusters with substructure scat- 
ter low on this plot. Some of these contain subclumps that 
are falling into the cluster for the first time and which 
have raised its mass without yet significantly altering the 
emission-weighted temperature; others are clusters which 
are undergoing a core-bounce after a merger and for which 
the temperature is temporarily slightly too low. 

It is apparent from the larger scatter in the upper panel 
of Fig. 1151 that the virial mass cannot be accurately deter- 
mined from a cluster's X-ray properties when spatial infor- 
mation is not available. However, we can reduce the scatter 
considerably by moving to a more compact radius (see Ta- 
ble 0. This has several advantages: the number of clusters 
showing substructure is reduced, the properties are mea- 
sured at radii that are more accessible to X-ray observations, 
and the substructure that is measured is more directly re- 
lated to subclumps that are interacting strongly with the 
intracluster medium and that will be influencing the X-ray 
properties. 

The lower panel in Fig. H5l shows the Tx-Msoo relation. 
As can be seen, the relation is much tighter than before, and 
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Figure 15. Cooling-flow corrected soft-band X-ray temperature 
versus mass. In the upper panel the outer radius of the cluster 
is -Rvii-; in the lower panel it is -Rsoo- The solid lines are fits to 
the data of a power- law relation T x oc M 60 . The circles repre- 
sent clusters for which the substructure statistic is greater than 
0.1; other clusters are indicated with triangles. Large symbols 
represent the resimulated clusters presented in this paper; small 
symbols are clusters from MTKP02. 



interestingly, the scatter of high substructure clusters from 
the relation is much-reduced. The circle that lies furthest 
above the solid line is Cluster 13 which has recently merged 
and whose temperature is still enhanced by the associated 
compression of the ICM. This cluster was even further from 
the line in the original plot (the topmost large triangle) but 
did not show up as having substructure because the sub- 
clump was too close to the centre. The circle that lies fur- 
thest below the line in the bottom panel is Cluster 5 which 
temporarily has a cooler temperature than the average due 
to core expansion following a recent merger. 

Note that, for most of the clusters, the mass decreases 
slightly in the move from i? v j r to R500, but for the two left- 
most clusters in the bottom panel (Clusters 7 and 14) it 
changes significantly. It can be seen from Table |5|that these 
are the two clusters with the highest degree of substructure 
at the end. What has happened is that while 7? v j r encloses 
the subclump causing the substructure 7?soo does not. The 
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Figure 16. Cooling-flow corrected bolometric X-ray luminosity 
estimated from emission in the soft band versus mass. In the 
upper panel the outer radius of the cluster is iJ v j r ; in the lower 
panel it is -Rsoo- The solid lines show power-law relations of the 
form Lx oc M 1,7 . The symbols have the same meaning as in 
Fig. El 



luminosity of these clusters also decreases significantly and 
they would be easily seen as bimodel in any X-ray observa- 
tion. For this reason, we omit them from the scatter statis- 
tics listed in Tabled] 

Having eliminated the principal cause of the scatter in 
the Tx-M relation, we again tested to see whether any of the 
properties considered earlier now correlate with the residual 
scatter, but with negative results. 

Fies ll6l and ll7l show equivalent plots for the luminosity- 
mass and luminosity-temperature relations. The former 
has a similar variation with substructure as seen in the 
temperature- mass relation. The latter shows no obvious cor- 
relation of scatter with substructure, consistent with our 
earlier observation that motion in the Lx-Tx plane during 
mergers tends to be parallel to the mean relation. However, 
in the upper panel of Fig. 1171 there is a concentration of 
clusters toward the bottom-left showing substructure. This 
reflects the fact that the sample was selected on the basis of 
virial mass and contains a number of clusters for which the 
merging is not complete. Conversely, in any X-ray luminosity 
or temperature-limited sample there will be a bias against 
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Figure 17. Cooling-flow corrected bolometric luminosity versus 
X-ray temperature estimated from emission in the soft band. In 
the upper panel the outer radius of the cluster is R v i r ; in the 
lower panel it is Rsoo- The solid lines show power-law relations 
of the form Lx oc Tx 2 ' 8 . The symbols have the same meaning as 
in Fig. EH 

high virial mass clusters with a high degree of substructure. 
This bias is reduced when properties are measured within 
-R500 instead. 



5 DISCUSSION 

In this study, the formation of 20 large clusters of galax- 
ies has been followed by resimulating the regions around 
them at high resolution. The X-ray properties were found 
to vary with time, driven mainly by the merging history of 
the clusters. Accretion during major mergers tends to in- 
crease a cluster's luminosity, pushing it up the Lx-Tx rela- 
tion, while between major mergers it slowly decreases, fol- 
lowing the movement of the Lx-Tx relation. In addition to 
this long-term variation, there are short-lived fluctuations 
associated with the merger itself. 

Over the period investigated clusters were said to be un- 
dergoing a major merger about 25-30 per cent of the time. 
However, to convert this to a number that can be compared 
to X-ray observations is very complex. Firstly, our defini- 
tion of a major merger was very ad-hoc and a more detailed 



comparision of possible definitions with X-ray observability 
would be required. In addition, there will be strong selection 
effects that favour high-redshift clusters undergoing tempo- 
rary boosts in X-ray luminosity. A detailed investigation of 
these biases will be required in order to correcctly interpret 
observations of the X-ray properties of clusters at high red- 
shift. 

As a subclump crosses R vir , the mass and the luminos- 
ity increase but the temperature stays roughly constant or 
decreases slightly (since the subclump will be cooler than 
the cluster). This causes the cluster to move below the Tx- 
M and Lx-M relations and so clusters with structure within 
R viT tend to be scattered low on the plots. For this reason, 
the mass within a smaller region such as -R500 correlates 
more closely with X-ray properties than does the virial mass. 

As the subclump moves through the ICM of the clus- 
ter, it compresses the gas in front of it. This causes adi- 
abatic heating and can be observed as a hot planar com- 
pression f ront. These feat u res h ave been observed, for ex- 
ample by iMaughan et aD J2003T) . Fig. 5 from their paper 
shows surface-brightness contours from a Chandra map of 
CI J0152.7-1357 a cluster at z = 0.833 and shows that the 
cluster is undergoing a major merger. Fig. 11 from the same 
paper shows the hardness of the X-rays with lighter regions 
representing harder radiation. Given that harder X-rays will 
be produced by higher temperature regions it can be seen 
that the gas between the merging clusters has been heated 
as it is compressed, similar to our clusters (see panel 1 of 
Figs. GO EH & ED- 

By the time the subclump reaches -R500 the cluster will 
be starting to undergo a merger boost and will start to move 
up in the Lx-M and Tx-M planes back toward the mean re- 
lations. This means that clusters with structure will not nec- 
essarily scatter off the Tx-M or Lx-M relations when prop- 
erties within this radius are considered. Usually the lumi- 
nosity and temperature increase together pushing the clus- 
ter roughly parallel to the Lx-Tx relation. However, if the 
merger triggers cooling within the cluster core then the X- 
ray temperature decreases and the cluster will scatter above 
the Lx-Tx relation and below the Tx-M relation. 

As the subclump moves toward the core the compres- 
sion front will get hotter and closer to the core of both the 
subclump and the cluster. This means that it will start to 
ram-pressure strip the subclump of its diffuse ICM. There- 
fore even if the subclump 's core survives the merger for a 
time, the diffuse gas will be assimilated into the cluster by 
the peak of the merger. 

Once the subclump has reached the core, the luminos- 
ity and temperature boost will have reached their maximum 
but, unless the infalling clump has a significant impact pa- 
rameter, the cluster wil l be roughly spheric al and appear to 
have little substructure. lGrego et all j2004) observe that the 
luminous, 2 = 0.783 cluster MS1137.5+6625 appears spher- 
ical in the optical and in X-rays but closer inspection with 
Chandra observations show that the cluster is not relaxed. 
The mass distribution is very compact, consistent with a 
large amount of recently accreted material. Since the lu- 
minosity can be increased at this time by up to an order 
of magnitude above what a truly relaxed cluster of similar- 
mass would be, then this could impose a bias in a flux lim- 
ited sample of clusters, particularly at high redshift, toward 
clusters at the point of merging. 
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The effect of merger boosts in biasing the estimate of 
cluster parameters from high-r edshift X-ray c l uster observa- 
tions has been investigated bv lRandall et all l|2002f ). 

After the peak of the collapse, the core of the sub- 
clump will continue to move past or through the core of 
the cluster. This behaviour has also been observed. Chan- 
dra observations of 1E 0657-56 (z — 0.296) presented by 
iMarkevitch et all i2002T) observe a 'bullet' which is surmized 
to be the core of a subcluster. The subcluster has previ- 
ously passed through the core of the main cluster (some 
0.1-0.2 Gyr ago) and has had its surrounding gas removed 
by ram-pressure stripping. It is unclear whether this object 
exceeds the escape speed of the cluster or not, but it still 
shows that a subclump can pass through the core of a cluster 
and that the haloes can merge while the two cores continue 
on their paths. 

If, as is likely, the core of the infalling subclump does 
not exceed the escape speed, then it will at some later point 
return to the core and merge. This will cause the lumi- 
nosity of the core to increase once more. If the tempera- 
ture profile does not change significantly when this hap- 
pens, then the increase in emission at the cooler core will 
cause an emission-weighted temperature for the whole clus- 
ter to decrease. This will move the cluster up and to the 
left on the Lx-Tx plane pushing the cluster above the Lx- 
Tx relation and so could contribute to the scatter about 
the relation. By this time the rest of the cluster will have 
settled down again. This will mean that the luminosity 
boost at the core will make the core much brighter than 
the surrounding temperature or luminosity should suggest. 
This could make the cluster look briefly like a cooling-flow 
cluster with a very high accretion rate. We find instan- 
taneous mass accretion rates M = L/ (5kT /2/imif) of up 
to 700 h~ 2 M@ yr -1 similar to o b served values in high red- 
shift clusters lEdee et al] 11994 iFabian fc Crawford! Il995t 
lAllen. Fabian fc Kneiblll996UAllen et aljll996ft ~ 

The current problem with cooling-flow clusters is that 
icooi <SC tn and s o a lot of gas should have cooled, but this 
is no t observed llKaastra et al.l 120011 : iPeterson et_alj j200lj. 
|2003|). Unfortunately, even considering the iLokenet^tH 
(1999) result that cooling-flow clusters usually occupy 
densely populated regions where clusters are more likely to 
be undergoing mergers, this explanation cannot explain the 
lack of cold gas in all cooling-flow clusters. A ROSAT sur- 
vey jPeres et aljll99sl) found that 70-90 per cent of clusters 
have cool cores. By their very nature the core mergers are 
short lived and so would be rare. However, it must be ac- 
cepted that the accretion of subclumps will continually feed 
the cores of clusters with low-entropy gas and this must be 
a contributory factor to the resolution of the cool-core prob- 
lem. The accretion of low-entropy material by clusters will 
be investigated in a future paper. 
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